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Background 

This invention generally relates to probabilistic models for calculating wire 
density in different areas of a datapath or hardmac, and more specifically relates to a 
probabilistic model which differentiates between horizontal and vertical segments. 
5 FIGURE 1 shows a general datapath structure. As shown, datapath cells are 

located in blocks (clusters) and control cells may be located in columns along left and 
i- right sides of the datapath and in blocks as well. Current datapaths are typically very 
- large and complex, and require design datapath macros with complex hierarchical 
= structures that have a high design quality. Additionally, complex constraints on the 
10 placement of cells, pins and nets, as well as the size of gaps between cells, blocks, etc. 
must be respected, and there must be guaranteed 100% detail routing. 

One prior art approach for calculating wire density in different areas of a 
datapath is based on a simplified density model and is used for placement quality 
estimation only. The approach is not accurate, does not differentiate between vertical 
15 and horizontal segments of connections, does not take into account all possible 

shortest length configurations of connections, and is unacceptable for calculating a 
congestion map. 
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Objects and Summary 

A general object of an embodiment of the present invention is to provide a 
method of estimating wire densities which differentiates between horizontal and 
5 vertical segments. 

Another object of an embodiment of the present invention is to provide a 
method of estimating wire densities which facilitates the formulation of a congestion 
- map. 

= Still another object of an embodiment of the present invention is to provide a 

10 probabilistic model which takes into account all possible shortest length 

configurations of connections, thereby being sufficiently accurate to estimate wire 
* density. 

Briefly, and in accordance with at least one of the foregoing objects, an 
embodiment of the present invention provides a method of accurately estimating 

15 horizontal and vertical wire densities in a datapath or hardmac. The method provides 
that the datapath or hardmac is divided into areas, and mathematical expectations are 
calculated for full and partial horizontal and vertical segments for each of the areas. 
The mathematical expectations are summed for both the horizontal and vertical 
segments, and this is done for each connection within the datapath or hardmac in order 

20 to estimate both horizontal and vertical wire densities. A congestion map can be 
created, and 100% detail routing is effectively guaranteed as a result of using the 
method. 



Brief Description of the Drawings 

The organization and manner of the structure and operation of the invention, 
together with further objects and advantages thereof, may best be understood by 
reference to the following description, taken in connection with the accompanying 
drawings, wherein like reference numerals identify like elements in which: 

FIGURE 1 is a general schematic diagram of a datapath structure; 

FIGURE 2 is a general schematic diagram showing a datapath area divided into 
square areas, and two pins, PI and P2, contained in the datapath area; 

FIGURE 3 is a general schematic diagram of a portion (i.e., rectangle [a,b,c,d]) 
of the area shown in FIGURE 2; 

FIGURES 4a-4f show six different possible connection configurations in a 
given area A; 

FIGURE 5 is similar to FIGURE 3, showing a full horizontal segment through 
area A in rectangle [a,b,c,d]; 

FIGURE 6 is similar to FIGURE 3, showing a partial horizontal segment 
through area A in rectangle [a,b,c,d]; 

FIGURE 7 is similar to FIGURE 3, showing another type of partial horizontal 
segment through area A in rectangle [a,b,c,d]; 

FIGURE 8 is a general schematic diagram showing a datapath area divided into 
square areas and a bounding area [a,b,c,d] defined in the datapath; 

FIGURES 9a and 9b show two different minimum bends configurations; 



FIGURE 10 shows bounding area [a,b,c,d] and the four different ways in which 
a minimal bends connection can go through a given area; 

FIGURES 1 la shows the horizontal probabilities for the configuration shown 
in FIGURE 9a; 

FIGURE 1 lb shows the horizontal probabilities for the configuration shown in 
FIGURE 9b; 

FIGURE 12 shows the overall horizontal probabilities for the configurations 
shown in FIGURES 9a and 9b; 

FIGURE 13 is similar to FIGURE 12, but shows the overall vertical 
probabilities; 

FIGURES 14a and 14b show two different two bends configurations; 
FIGURES 15 shows the horizontal probabilities for the configuration shown in 
FIGURE 14a; 

FIGURE 16 shows the horizontal probabilities for the configuration shown in 
FIGURE 14b; 

FIGURE 17 shows the overall horizontal probabilities for the configurations 
shown in FIGURES 14a and 14b; and 

FIGURES 18a and 18b show two different three bends configurations. 
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Description 

While the invention may be susceptible to embodiment in different forms, there 
are shown in the drawings, and herein will be described in detail, specific 
embodiments with the understanding that the present disclosure is to be considered an 
5 exemplification of the principles of the invention, and is not intended to limit the 
invention to that as illustrated and described herein. 
9 An embodiment of the present invention provides a probabilistic method for 

1 calculating wire density in different areas of the datapath (the term "datapath" is to be 
% construed very broadly herein, as the term is used herein to mean any type of real 
-10 estate) and other hardmacs with a given cell placement. The method is based on a 

probabilistic model of connection between two pins. The model takes into account all 
C possible shortest length configurations of the connection, and differentiates between 
vertical and horizontal segments of the connection. Thus, the model is sufficiently 
accurate to be used for wire density estimation, and provides that congestion maps can 
15 be calculated. 

Initially, as shown in FIGURE 2, a datapath area is divided into M DP by N DP 
squared areas (effectively a matrix of columns and rows), where each area size is 
about equal to the width of placement columns or cell width. It is known that the 
number of the shortest length paths (configurations) from PI to P2 (see FIGURE 3) is: 

20 

<» ■G,°-i a )- (m ( °y-^ l 1) , 
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With, regard to the possible shortest length configurations for connection C 
from pin PI to pin P2 (see FIGURE 1), it is known that any shortest length connection 
(PI, P2) will belong to rectangle [a,b,c,d]. Therefore, for any area Al from [a,b,c,d] 
there exists some probability P(A1)>0 that connection (PI, P2) will go through this 
area, and for any area A2 outside [a,b,c,d] there is zero probability that connection (PI, 
P2) will go through this area. If the probability P(a) is known, then the mathematical 
expectation of area A having connection (PI, P2) can be calculated as being M(A) = 
P(A). 

As shown in FIGURES 4a-4f, any connection (PI, P2) can go through area A in 
six different ways. The configurations shown in FIGURES 4a-4d are possible if pin 
P2 is higher than pin PI, and the configurations shown in FIGURES 4e and 4f are 
possible if pin PI is higher than pin P2. The segment (i.e., the connection through 
area A) shown in FIGURE 4a is a full horizontal segment, while the segment shown in 
FIGURE 4b is a full vertical segment. Each one of the segments shown in FIGURES 
4c-4e is both a partial horizontal and partial vertical segment. 

The situation where pin P2 is higher than pin PI (as shown in FIGURE 2) will 
be discussed below. The analysis would be analogous in the case where pin PI is 
higher than P2. 



FIGURE 5 shows the rectangle [a,b,c,d] (see FIGURE 2) wherein there is a full 
horizontal segment (see FIGURE 4a) through area A (FIGURES 5-7 also show the 
numeration of the columns and rows which form the rectangle [a,b,c,d])- The 
mathematical expectations M hl (A) of full horizontal segments can be calculated as 
5 follows: 



(2) M (A) - N(P1,A').N(A",P2) 

{ ) Mhl(A) ~ N(P1,P2) 



where 

N(P1, A') is the number of possible paths from PI to area A'; 
N(A", P2) is the number of possible paths from area A" to P2; and 
N(P1, P2) is the number of possible paths from PI to P2. 



Taking into account the formula (equation (1) above) for the number N(P1, P2) 
of paths between two points, all the numbers which are needed to calculate the 
mathematical expectation M hl (A) can be found: 
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(4) 



FIGURE 6 shows the rectangle [a,b,c,d] (see FIGURE 2) wherein there is a half 
or partial horizontal segment (see FIGURE 4c) through area A. The mathematical 
expectations M h2 (A) of half horizontal segments can be calculated as follows: 



f*\ M N(P1,A')-N(A B, ,P2) 
(5) M h2 (A) = 0.5 N(P1,P2) 



where 

Z N(A'", P2) is the number of possible paths from area A'" to P2, and coefficient 0.5 
1 0 indicates that there is only half of a horizontal segment in area A. 

Taking into account the formula (equation (1) above) for the number N(P 1 , P2) 
* of paths between two points, all the numbers which are needed to calculate the 

mathematical expectation M h2 (A) can be found: 

15 (6) N(A'", P2> = (m-i^-j-l)t 

(m-i-1)! - (n-j)! 
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FIGURE 7 shows the rectangle [a,b,c,d] (see FIGURE 2) wherein there is 
another type of half or partial horizontal segment (see FIGURE 4d) through area A. 
The mathematical expectations M h3 (A) of half horizontal segments can be calculated 
as follows: 



where 

10 N(P1, A"") is the number of possible paths from PI to area A"", and coefficient 0.5 
indicates that there is only half of a horizontal segment in area A. 

Taking into account the formula (equation (1) above) for the number N(P1, P2) 
of paths between two points, all the numbers which are needed to calculate the 
1 5 mathematical expectation M h3 (A) can be found: 
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(7) 



M h3 (A) = 0.5 - 



NfPl,A"")-N(A,P2) 
N(P1,P2) 



(8) 



N(P1, A"") 



(i-2)!-(j-l)! 
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To determine the whole mathematical expectation M h A1, (A) of all horizontal 



segments of all connections, the following summation is calculated: 



(9) 



c e Connections 
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where M h c (A) = M h (a) is the whole mathematical expectation of horizontal segments 
- in area A for one connection c. 

The same approach can be used to obtain formulas for vertical segments and 
10 the case where PI is higher than P2. 

From the formulas above, it can be concluded that the time complexity of the 
algorithm will depend on how fast factorials (n!) Can be calculated. If a 
straightforward calculation is used, then the time complexity for one connection and 
15 one area is 0(m+n). The time complexity for one connection and all areas (see 

FIGURE 2) is 0(M DP N DP (m+n)). Finally, the time complexity for all N connections 
and all areas is 0(M DP N DP N(m+n)). There are several ways how to deduce the time 
complexity, especially for long connections with large m and n. One method is to use 
the Sterling formula for factorial calculation: 
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(10) 



n! 
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Then, the time complexity becomes 0(M DP N DP N). The same time complexity and 
even better time can be obtained if factorials of integer numbers are tabulated in 
advance for the range of approximately [1-100]. 

The time efficient method (time complexity is proportional to the product of the 
connections and areas) described above can be used to accurately estimate horizontal 
and vertical wire density in different areas of datapath or hardmac. The approach is a 
good probabilistic model for connections going through areas with high wire density. 
The model differentiates between horizontal and vertical segments, and takes into 
account all possible shortest length configurations of connections. The model also 
provides for the calculation of a congestion map. 

However, the approach described above has the following two drawbacks for 
chip areas with low and middle wire density. First, it assumes that the connection can 
have any configuration with the same probability. This is not always true as the 
connection more likely has a configuration with a small number of bends in chip areas 
with low and middle wire density. Second, it assumes that the probability of any 
connection configuration that goes through or near the center of the bounding box (i.e. 
rectangle [a,b,c,d]) (see FIGURE 8) around the connection is higher than for other 
configurations. This is not always true as it will depend on the location of other pins 
and wires. 
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A better approach uses the model with minimum bends in areas with low wire 
density, and uses models with more bends in areas with middle and high wire density. 
The rule is: "the more wire density the more bends in the model". First, the model 
with minimum bends is found, then the model is used recursively to build other 
models with more bends. 

Initially, the chip is divided into M DP by N DP squared areas as shown in 
FIGURE 8, where each area size is about equal to the width of placement columns (or 
cell width) (in some cases the chip may be rotated 90 degrees, hence there are no 
placement columns, but rather placement rows). 

The minimum bends model (i.e. model 1) describes all connection 
configurations with only one bend and the shortest length. FIGURES 9a and 9b show 
these configurations for connection C from pin PI to pin P2. 

The probability P(a) for each area A of the connection bounding box [a,b,c,d] 
to have the connection (PI, P2) will now be found (see FIGURE 10). For any area Al 
from [a,b,c,d] there exists some probability P(A)>0 that connection (PI, P2) will go 
through this area, and for any area A outside [a,b,c,d] there is zero probability that 
connection (PI, P2) will go through this area. If probability P(A) is known, then the 
mathematical expectation of area A having connection (PI, P2) is M(A) = P(A). Any 
connection (PI, P2), where PI is lower than P2, con go through area A in four 
different ways as shown in FIGURE 10. The case where pins PI and P2 are placed as 
shown in FIGURES 9a and 9b will be considered. For the situation where pin PI is 
higher than pin P2, the analysis will be analogous. 
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FIGURE 1 la shows the bounding box [a,b,c,d] (see FIGURE 10) (FIGURES 
10, 1 la, 1 lb, 12, 13 and 15-17 also show the numeration of the columns and rows 
which form the rectangle or bounding box [a,b,c,d]), and all horizontal probabilities 
(mathematical expectations) for the configuration shown in FIGURE 9a. Areas with a 
full horizontal segment have 0.5 probability, because there are only two possible 
configuration, while areas with a half horizontal segment have 0.25 probability, 
because these areas contain about 0.5 part of the segment and there are only two 
possible configurations. FIGURE 1 lb shows all the horizontal probabilities for the 
configuration shown in FIGURE 9b. 

The whole mathematical expectation M h (A) can be found as a sum: 

(11) M h (A) = M hl (A) + M h2 (A) 

of mathematical expectations for both configurations shown in FIGURE 1 1 (see 
FIGURE 12). The formula for mathematical expectation M h (A) is as follows: 

(12) M h (A) = 0.5 if i = 1 andj = 2, 3,...,n-l 

(13) M h (A) = 0.5 ifi = mandj = 2,3,...,n-l 

(14) M h (A) = 0.25 if i= 1 and j = 1 or j = n 

(15) M h (A) = 0 if i=2, 3,...,m-l and j = 1, 2,..., n, 
where local (inside [a,b,c,d]) numeration of rows and columns is used. 
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The same formulas can be used for horizontal segments when point PI is higher 
than point P2. 

To determine the mathematical expectation M h AI1 (A) of all horizontal segments 
of all the connections, the following summation is calculated: 

(16) Mf(A)= J M c h (A) 

c € Connections 

where M h c (A) = M h (A) is the whole mathematical expectation of horizontal segments 
in area A for one connection c. 

The same approach is used to obtain formulas for vertical segments (see 
FIGURE 13): 

(17) M V (A) = 0.5 ifj = landi = 2,3,...,m-l 

(18) M V (A) = 0.5 ifj=mandi = 2,3,...,m-l 

(19) M V (A) = 0.25 if j = 1 andi=l ori = m 

(20) M V (A) = 0 if j = 2, 3,..., n-1 and i - 1, 2,..., m, 

where local (i.e. inside [a,b,c,d]) numeration of rows and columns is used. 

The same formulas can be used for vertical segments when point PI is higher 
than point P2. 
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To calculate the mathematical expectation M V AI! (A) of all vertical segments of 
all the connections, the following summation is calculated: 

(21) Mf (A) = J Ml(A) 

c e Connections 

5 

From the formulas above, it can be concluded that the time complexity of the 
model will depend on n and m. The time complexity for one connection is 0(m+n). 
== The time complexity for all N connections is 0(N(m+n)). 

Next, the obtained formulas for one bend configurations are recursively used to 
10 find models with 2, 3 ... bends. With regard to a connection configuration with two 
r bends, there are two possible types of configurations, and these are shown in 

FIGURES 14a and 14b. 

To determine all mathematical expectations, two bend configurations are 
considered as a combination of all possible one bend configurations (PI, P2'). There 
15 are m possible locations for P2' for the configurations shown in FIGURE 14b. The 
whole mathematical expectation M h (A) can be found as a sum: 

(22) M h (A) = M hl (A) + M h2 (A) + . . . + M h(m+n) (A) 

20 of mathematical expectations for all possible configurations in FIGURE 14. FIGURE 
15 shows the mathematical expectations for the configuration shown in FIGURE 14a, 
while FIGURE 1 6 shows the mathematical expectations for the configuration shown 
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in FIGURE 14b. FIGURE 17 shows the overall horizontal mathematical expectations. 



The formula for mathematical expectation M h (A) is as follows: 



(23) 


M h (A) 


= 0 . 5(mn+n-(j- 1 )m)/nm 


ifi 


= 1 andj =2, 3,..., n-1 


(24) 


M h (A) 


= 0.5(jm+n)/nm 


ifi 


= mandj =2, 3,..., n-1 


(25) 


M h (A) 


= 0.25(m+l)/m 


ifi 


= 1 and j = 1 


(26) 


M h (A) 


= 0.25(n+m)/nm 


ifi 


= m and j = 1 or j = n 


(27) 


M h (A) 


= 0.25(m+l)/nm 


ifi 


= m and j = n 


(28) 


M h (A) 


= 0.5/m 


ifi 


= 2, 3,..., m-1 and j = 1, 2,..., n, 



where local (i.e. inside [a,b,c,d]) numeration of rows and columns is used. 



The same formulas can be used for horizontal segments when point PI is higher 
than point P2. The same approach can be used to obtain formulas for vertical 
segments: 



(29) 


M V (A) = 0.5(mn+m-(i-l)n)/nm 


ifj 


= 1 and i 


= 2, 3,..., m-1 


(30) 


M V (A) = 0.5(in+m)/nm 


ifj 


= n and i 


= 2, 3,..., m-1 


(31) 


M V (A) = 0.25(n+l)/n 


ifj 


= 1 and i 


= 1 


(32) 


M V (A) = 0.25(n+m)/nm 


ifj 


= n and i 


= 1 or i = m 


(33) 


M V (A) = 0.25(n+l)/nm 


ifj 


= n and i 


= m 


(34) 


M V (A) = 0.5/n 


ifj 


= 2,3,..., 


n-1 and i = 1, 2,..., m, 



where local (i.e. inside [a,b,c,d]) numeration of rows and columns is used. 
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From the formulas above, it can be concluded that the time complexity of the 
model will depend on n and m. The time complexity for one connection and one area 
is O(mn). The time complexity for all N connections is 0(Nmn). With the increase of 
bends in the model, the time complexity also increases. 

5 A three bends model will now be outlined. To arrive at the three bend model, 

the obtained formulas for the two bend configuratiosn will be used. FIGURES 1 8a 

I and 1 8b show the two possible three bend configurations. To determine all the 

mathematical expectations, the three bend configurations are considered as being a 
combination of all possible two bend configurations (PI, P2'). There are n possible 
TO locations for P2' for the configuration shown in FIGURE 1 8a, and m possible 

locations for P2' for the configuration shown in FIGURE 1 8b. Using the approach 
described above recursively, a model with any given number of bends can theoretically 
be built. However, the calculations to build k-bends (k>3) models may be such that 
such a model would be impractical in light of the large expense compared to the 
15 relatively small improvement in accuracy. To increase the speed of all the 

calculations, all possible matrices for the mathematical expectations for given sizes of 
connections m and n can be tabulated (there will be mn matrices for each model). 
Then, for any k-bends model, the time complexity for all N connections will always be 
O(Nmn) due to the fact that the time for tabulation can be ignored since the tabulation 
20 is only perormed once. 
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The time efficient models described above can be used to accurately estimate 
horizontal and vertical wire density in different areas of datapath or hardmac. These 
models take into account all possible minimal bends and shortest length configurations 
of the connection. Thus, these models are accurate enough to be used for wire density 
estimation in areas with low, middle and high wire density. Preferably, the model with 
minimum bends is used in areas with low wire density, and models with more bends 
are used in areas with middle and high wire density. 

While embodiments of the present invention are shown and described, it is 
envisioned that those skilled in the art may devise various modifications of the present 
invention without departing from the spirit and scope of the appended claims. 
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